Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN M?eX style file vl.4) 



Hot dust in two hard Chandra X-ray sources 



R.J. Wilman, A.C. Fabian and P. Gandhi 

Institute of Astronomy, Madinglcy Road, Cambridge CBS OHA 



1 February 2008 



o 
o 
o 

< 



OO 
00 

o 
o 
o 

6 



ABSTRACT 

The two brightest hard X-ray sources discovered serendipitously by Chandra in the 
field of the lensing cluster A2390 are found to have /S'OCAM counterparts at 6.7 and 
15 /im. We use this fact, together with their non-detection by SCUBA at 850 /im, as 
the basis for DUSTY radiative transfer modelling of their infrared spectral energy 
distributions. For the best-fitting models, we find that the dust which reprocesses 
the optical-UV light in these Compton-thin AGN is heated to near its sublimation 
temperature (above 1000 K), with an inner radius within a parsec of the nucleus. Some 
warm dust models with inner temperatures of 200 K are also formally acceptable. 
These findings strongly support the obscured AGN hypothesis for the new hard X-ray 
Chandra sources, which lack both strong emission lines and 850 /im SCUBA detections. 



Key words: galaxies: active - quasars: general 
X-rays:general 



galaxies :Seyfert - infrared: galaxies ■ 



1 INTRODUCTION 

The Chandra X-ray observatory has recently resolved the 
bulk of the 2-7 keV hard X-ray background (HXB) into 
point sources (Mushotzky et al. 2000; Brandt et al. 2000). 
Much effort is now focused on determining the physical 
nature of these new sources and their relation to exist- 
ing classes of active galactic nuclei (AGN). Mushotzky et 
al. (2000) found that roughly a third of their sources were 
blue broad- line (type I) AGN, one third were identified with 
apparently normal galaxies, whilst the final third had an 
extremely faint or no optical counterpart (/ ^ 23mag) (see 
also Fiore et al. 2000). Crawford et al. (2000) confirmed 
that sources in the last group are easily detected in the 
near-infrared, with featureless spectra suggesting that the 
AGN is heavily obscured and colours consistent with red- 
dened elliptical galaxies at z — 1 — 2. Further insight has 
come from cross-correlating Chandra source lists with those 
from 850 fim SCUBA observations of the same fields: for the 
lensing clusters A2390 and A1835, Fabian et al. (2000) (here- 
after FOO) found only one source common to both datasets, 
implying that any AGN in the SCUBA sources must either 
contribute little to the sub-mm power or be Compton-thick 
with X-ray scattering fractions below 1 per cent; this result 
was confirmed by Hornschemeier et al. (2000), who detected 
none of the 10 sub-mm sources in and around the Hubble 
Deep Field North (HDFN) in a 166-ks Chandra exposure. 

The two brightest X-ray sources in the A2390 field, CX- 
OUJ215334.0+174240 and CXOUJ215333.2-H 74209 (here- 
after sources A and B respectively), have /SO CAM counter- 
parts at 6.7 and 15 ^m in the literature (Altieri et al. 1999; 
Lemonon et al. 1998). Together with their non-detection at 



850 /^m by SCUBA, this has implications for the location 
and properties of the dust which may be associated with the 
X-ray obscuration. In this paper we use the dust radiative 
transfer code DUSTY (Ivezic, Nenkova & Elitzur 1999) to 
model the infrared spectral energy distributions (SEDs) of 
these two sources as reprocessed optical-UV emission from 
an AGN. The faint optical light from these sources {V ~ 25) 
is assumed to be from the host galaxy and the HYPERZ 
code (Bolzonella, Pello & Miralles 2000) is used to estimate 
photometric redshifts. 



2 SED MODELLING 

For the primary optical-UV AGN continuum input to the 
DUSTY models, we follow Granato, Danese & Frances- 
chini (1996) and use a broken power-law with a = —0.5 
for log u < 15.4, Q = -1.0 for 15.4 < log u < 16, and 
a — —2.2 for log u > 16 (L^ oc The central source is 
surrounded by a spherical dust cloud, specified by its optical 
depth T at 0.3 fim {Ay ~ 0.57t, and for a Galactic dust:gas 
ratio A^H — 1.1 x lO^^rcm"'^), the ratio R = rout/nn of 
the radii of the outer and inner edges of the distribution 
(between which we take the density to be constant), and 
the dust temperature Tin at rin. The code currently only 
supports grains of a single type, so we adopt the default 
'standard ISM mixture' comprising grains which mimic a 
53:47 mixture of silicate and graphite grains from Draine 
& Lee (1984). The DUSTY wavelength grid spans 0.01 jum- 
3.6 X 10* fj,m, so the effects of harder radiation are not mod- 
elled. Since sources A and B are likely to be Compton-thin 
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ratios in FOO and from AGN synthesis models for the XRB, 
e.g. Wilman & Fabian 1999), they are transparent to pho- 
tons above a few keV; there are thus about 1.5 decades of 
frequency below 0.01 /im where photons would be absorbed 
but which arc not included. Given that the input SED cuts 
off sharply below 0.03 jj,m, this omission is unimportant. 

We adopt the cosmological parameters Ho = 
50kms~^ Mpc"'^ and qo = 0.5 throughout the paper. 

2.1 Source B (CXOUJ215333.2+174209) 

Lemonon et al. (1998) give 6.7 and 15 /xm 750 CAM fluxes 
for source B of llOl^gQ and 350l4[] //Jy, respectively. From 
a total of 33 counts, FOO estimate 0.5-2 and 2-7 keV fluxes 
of 5.9 and 23xl0~^^ ergcm"^ s~^ , respectively; the upper 
limit at 850 /xm is 5.7 mjy (FOO). We use the relative fluxes 
in the B, V(F555W), R, I(F814W), J and K' filters provided 
by Lemonon et al. along with HYPERZ to compute a pho- 
tometric redshift (N.B. the HST F555W and F814W magni- 
tudes given by Lemonon et al. diflfer from those listed in FOO, 
but for internal consistency we adopt the former values) . We 
consider 2 families of HYPERZ models: (i) where at least 80 
per cent of the K' light (2.103 /xm) is host galaxy starlight; 
(ii) where at least 50 per cent of it is nuclear dust emis- 
sion. A satisfactory HYPERZ model could not be found for 
the case whore all the K' light is from the nucleus. HYPERZ 
uses Bruzual & Chariot (1993) spectral synthesis models and 
varies the redshift, age and extinction of the population. The 
elliptical galaxy model we use has an exponentially-declining 
star formation rate with an e-folding time of 1 Gyr; other 
models with longer star-formation timoscalcs (SO-Sd and 
Im), as well as bursts, were also considered, but they do not 
significantly change the fitted redshift (essentially because 
the break between the V(F555W) and R points is assumed 
to be the redshifted 4000A break). 

For case (i), HYPERZ fits a 500 Myr old elliptical 
galaxy with Ay = 2.40 mag at z = 0.575. For this red- 
shift, we generated a grid of DUSTY models, with r = 
10, 20,. .,100 (equivalent to JVh ~ 10^^ - 10^^ cm'^ for 
a Galactic dustigas ratio, and appropriate for the mod- 
erately Compton-thin obscuration implied by the X-ray 
fluxes), Tin = 100, 200,. ., 1500 K, and R = 5, 50, 250 and 
1000. After integration over the ISOCAM filter bandpasses, 
the models were normalised to the observed 15 jira flux 
density, and deemed acceptable if the 6.7/15 ^ra flux ra- 
tio fell within the range allowed by the errors on the 
ISOCAM data, if the 850 /im flux feU below the SCUBA 
limit, and if the K' flux was less than 20 per cent of 
that observed. Fig. 1 shows which models moot those cri- 
teria, and Fig. 2 the SED fits for three such cases, with 
(Ji,T,Tin) =(250,50,1500), (50,100,700) and (5,60,200). The 
implied optical-UV DUSTY input luminosity, 2-10 keV X- 
ray absorption corrected luminosity (derived from the ob- 
served 2-7 koV flux, assuming an intrinsic power-law spec- 
trum with a photon index F = 2), optical (2500A) to X-ray 
(2keV) spectral index, Oox, and the inner radius, rin, for 
each are shown in Table 1. The best fits to the portion of 
the SED covered by the ISOCAM data are obtained with 
hot dust (heated close to its sublimation temperature of 
1500 K), with an inner radius within a parsec of the central 
engine; the implied Qox in this case is, however, somewhat 
flatter than the canonical value of 1.35 for quasars (Elvis et 
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Figure 1. The points show which of the DUSTY grid models 
described in section 2.1 for source B case (i) are acceptable, 
for different values of R = rout/nn. 

al. 1994). Figs. 1 and 2 also show, however, that models with 
warm dust at Tin = 200 K cannot be ruled out at present; 
such models have inner radii of tens of pc and Qox values 
close to those of quasars. Discrimination between the warm 
and hot dust models would be possible using data at 70 /im; 
e.g. the 200 K model shown in Fig. 2 has a 70 /im flux den- 
sity of 16 mJy, which is well within the capability of SIRTF 
(see e.g. Brandl et al. 2000). 

For case (ii), HYPERZ fits a 3.5 Gyr old elliptical galaxy 
at z=0.505 with Ay = 0.60 mag. The need for DUSTY to 
reproduce at least half of the K' light in this case moans 
that the dust must be hot (Ti„ = 1000 - 1500 K) and not 
too optically thick at 2 /xm (r ~ 25 — 40); R is not very well 
constrained, with a value of 100 for r — 35, although a more 
compact structure, (R,t) — (5,40), is also acceptable. For 
{R,t) = (100,35), we find L2-10 = 4.6 x lO-^^ergs-^ and 
aox = 1.12; the latter model SED is shown in Fig. 2. 

We conclude that the best-fitting models are those in 
which the dust is hot, with an inner radius within 1 pc of 
the nucleus. Several models with warm dust (200 K) arc also 
formally acceptable, but for R = 50 some of these underpre- 
dict the SCUBA upper limit by only a small factor (~ 2), 
which, if this source is representative of the new Chandra 
HXB sources, seems unlikely given the non-detection of large 
samples of them by SCUBA (see references in section 1). 

2.2 Source A (CXOUJ215334.0-|-174240) 

This is the brightest Chandra source in the A2390 field, co- 
incident with a slightly reddened L* mid-type spiral at a 
photometric redshift of 2; = 0.85 ± 0.15 (FOO). From a total 
of ~ 90 counts, FOO performed crude X-ray spectral fitting 
and deduced an intrinsic iVu = (6±2) - (9±3) x 10^^ cm"^ , 
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Table 1. Details of three acceptable DUSTY models for source B case (i), spanning the range of allowed temperatures. 





Model 1 


Moaei 2 


Moaei 3 


T 


50 


100 


60 


Tin (K) 


1500 


700 


200 


R = r-out/rin 


250 


50 


5 


nn (pc) 


0.41 


3.9 


67 


L{opt-UV) (1044 erg s-1) f 


6.2 


11.3 


19.2 


L(2 - 10 keV) (1043 ergs-i ) J 


6.3 


7.5 


6.5 




1.17 


1.24 


1.35 



f Luminosity of the DUSTY input continuum. 
I Absorption-corrected 2-10 keV luminosity. 
* Implied optical to X-ray spectral index of the intrinsic continuum. 
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Figure 2. Observed and model SEDs for source B, for the 
cases where the galaxy produces essentially all (case i) and 
half (case ii) of the light at 2.1 /im. The circles and trian- 
gles are the optical and ISOCAM fluxes, respectively, from 
Lemonon et al. (1998); the arrow is the SCUBA upper limit 
from FOO. The DUSTY models are normalised to the ob- 
served 15 nm flux density; the optical points are overlaid 
with the best-fit HYPERZ model for the host galaxy. In 
the upper panel, the solid, dashed and dotted lines are 
the models with (i?,r,Tin) =(250,50,1500), (50,100,700) and 
(5,60,200), respectively. 



with a de-absorbed L2-10 ~ 2 — 3 x 10'*'* ergs~* (after cor- 
rection for a lensing magnification by a factor of about 2), 
for z — 0.7 — 1, respectively. It is thus an X-ray type II 

quasar. 

It also has /SOC^M counterparts at 6.7 and 15 /im (Al- 
tieri et al. 1999) but the source fluxes have not been pub- 
lished. We proceed on the assumption that they are equal to 
those of source B, which appears reasonable from the images 
in Altieri et al. 

For redshifts in the range 0.7-1.0, we find acceptable 
DUSTY models with r = 55 - 80 (equivalent to the fitted 
A'h = 6 — 9 X 10^^ cm~^ for a Galactic dustigas ratio) for 
R = 100 and Tin = 1500 K; they have nn = 0.55-1.0 pc and 
Qox ^ 1.0. 



3 DISCUSSION 

It is instructive to compare our results with the calculations 
of Granato, Danese & Franceschini (1997), who compared 
dust radiative transfer models with the infrared properties of 
Seyfert galaxies. They concluded that the moderately thick, 
extended tori of Granato & Danese (1994) (GD) (with 5 < 

< 80 mag and outer radii of lO's to lOO's of parsec), 
provided a better fit to the data than the thick, very compact 
models of Pier & Krolik (1992a) (PK) (with Ay > 800 mag, 
all within a pc). They also found that the observed X-ray 
absorption to the nucleus was much higher than that implied 
by the Ay obtained by modelling the infrared SED (for a 
Galactic dustigas ratio), suggesting that much of the X-ray 
absorbing gas lies within the dust sublimation radius. 

In terms of its optical depth and radial extent, the 
obscuring dust in these two Compton-thin, absorbed, X- 
ray background sources more closely resembles the model 
of GD than that of PK. Indeed, Granato, Danese & 
Franceshini (1997) predicted a significiant correlation be- 
tween the HXB sources and those appearing in ISO surveys 
at 10-20 /nm, precisely because the GD models are (at the 
most) only moderately optically thick over this wavelength 
range. An important difference, however, is that the GD tori 
have opening angles of 35-45 degrees, whereas the covering 
fraction of any torus-like structure in the HXB sources must 
be higher (e.g. Fabian & Iwaweisa 1999 deduce that 85 per 
cent of the accretion power must be absorbed). Note also 
that the DUSTY calculations assumed a spherical dust ge- 
ometry. It is thus appropriate to examine how such a high 
space covering obscuration could be maintained against the 
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dissipative forces which would tend to yield a flattened struc- 
ture: at the sub-parsec scale inner radii found for sources A 
and B, Pier & Krolik (1992b) have shown that the pres- 
sure of the nuclear radiation on the grains can make the 
torus geometrically thick, by both reducing the vertical com- 
ponent of gravity, and if the torus is sufficiently clumpy, 
driving random motions of its constituent clouds. Fabian 
et al. (1998) demonstrated that a space-covering obscura- 
tion could be maintained with energy input from a nuclear 
starburst within 100 pc, but the distinct lack of emission 
lines in manv of the newlv-discovered HXB sources appears 



to rule [this out. More generally, the lack of AGN emission 
lines in such objects may be due to dust within the narrow 
line region clouds. Netzer & Laor (1993) showed that its 
presence could significantly suppress line emission, through 
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tween the inferred covering factors of the broad and narrow 
line regions in classical AGN when dust is not considered 
(the broad line region lies just within the dust sublimation 
radius) . Our findings suggest the presence of a high covering 
fraction of dusty gas in the new HXB sources. 

Alternatively, the lack of emission lines and non- 
detections by SCUBA could be used to argue that the new 
HXB sources are not actually obscured AGN, but some other 
class of intrinsically hard source (e.g. ADAFs, as proposed 
for the HXB by Di Matteo & Allen 1999); the present quality 
of the X-ray spectra is not high enough to discriminate be- 
tween these two possibilities. Our present findings, however, 
strongly support the obscured AGN hypothesis, by demon- 
strating that the energy which is inferred to be absorbed 
in the optical-UV-X-ray range is reradiated in the infrared. 
In consequence, the results of deep Chandra surveys of well- 
studied ISO fields (e.g. from the FLATS consortium; Rowan- 
Robinson et al. 1999), are eagerly awaited. Indeed, as noted 
by FOO, 5 of the 6 Chandra sources discovered by Horn- 
schemeier et al. (2000) in the HDFN also have ISOCAM 
counterparts at 15 /xm listed by Aussel et al. (1999) (CXO- 
HDFN 123648.1+621309 is the only source not detected). 
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